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PREFACE 


During the final period of this grant a theoretical and experimental 
study was made of the effect of diffusion on the saturation of the infrared 
absorption of radiation by gases. This report contains most of the 
theoretical development of a model of a gas absorber including diffusion 
and a small amount of the experimental evidence we gathered to support 
this model. 

The experimental results along with the theoretical model will be 
submitted for publication in the near future. 



Abstrac t 


Effect of diffusion on the saturation and absorp- 
tion of the molecular gases is considered.. The analysis 
starts with a simple derivation of the multi-level rate 
equations. Inco rporating the diffusion term, calculations 
are carried out both for the homogeneous and inhomogeneous 
broadening cases. Both graphical and analytical solutions 
are obtained for the absorption coefficient and the satu- 
ration intensity. It is shown that the effect of diffu- 
sion is characterized by a diffusion function 
which varies inversely with the radius of the beam at low 
pressures and small beam sizes.. The analysis is done in 
terms of two parameters k]_a and kga, which are simple 
functions of the intensity of the beam I, the effective 
relaxation rate w, and the effective diffusion constant 
D. M a'* is the radius of the beam. 
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I. Introduction 


Absorption and saturation of infrared radiation 

around 10. 6ju in molecular gases has recently been of 

1 , 2 , 3,k 

considerable interest. Using a two-level model 

where one is only concerned with the lower and upper 
levels used in absorption, the variation of the laser 
beam intensity, I, is found to be 
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dl = - I (1) 

dz 

where ©£ is the absorption coefficient which can be ex- 
pressed as 

<=*= °*o (2) 

fr-fc-r 


In this expression I s is the saturation intensity; 
od c is the small signal absorption coefficient; n is 1 
when the transition is homogeneously broadened and . £ 
when the transition is inhomo gen eou sly broadened. 

If we define 

M = + M 2 

AM = M 2 - M x (3) 

m = AM/M 

where M^ and M 2 are the lower and upper level population 
equation (2) can be written as^ 
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1 + 



(k) 


Even though two-level model gives a very simple 
explanation of saturation, it is not totally satisfactory 
because of various phenomena that are involved in satu- 
rable gases. The most important effects are the tight 
coupling between various vibrational-rotational levels 
and the spatial diffusion of molecules. If multi- 
rotational levels are considered, C.P. Christensen et al*^ 
have showed that the results of the two-level model can 
be used by replacing relaxation rates and diffusion con- 
stants by their effective values as follows: 
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In these expressions w-^ and w ^ are the lower and 
upper effective relaxation rates; the ith rota- 

tional level relaxation rate in the upper(lower) band; 
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u . 1 

m^ is the corresponding population density; M-^ and 
are the lower and upper absorption level population den- 
sities. Similarly, and are the effective diffusion 
constants, and D is the diffusion constant. Assuming the 
rotational constant B for the upper and lower bands are 
the same, and will be nearly equal. A more simple 
derivation of these expressions than was done originally 
is included in the appendix. 

With the diffusion term included, the rate equation 
for the inversion density can be obtained from equation 
(4) as 6 

DV^n - (2o-I + w)m = w (6) 

where D is assumed to be the effective relaxation rate 
for both bands; w is equal to WgJ or" is the absorption 
cross-section at the laser wavelength, and I is the 
intensity of the beam in photons/cm^-sec. Since the 
lower level for most saturable absorbers is the ground 
state, w-^ is assumed to be zero. 

The absorption coefficient can be written as 

<=*=-1 dl = m <=* ( 7 ) 

—y — o 

I dz 

where i s the small signal absorption coefficient. 

Since all the variables are a function of the radial 
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coordinate r, we would like to define average quantities. 
First of all, average intensity I_„ is defined as 

v Q. V 

P = 7ra^ I (8) 

av 

where P is the power and a is the radius of the beam. 

Then equation ( 7 ) can be transformed to 


av 


- JL 

I 


dl 


av 


av 


dz 


= °<o m av 


where 
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av 


oo 

-'So 


-r^/ , 
me 2TTrdr 


( 9 ) 


( 10 ) 


for a Gaussian beam. If the beam is rectangular instead 
of Gaussian as discussed in the next section, the equa- 
tions would be the same except for 


m. 


av 


S a- 

0 " m 2 IT rdr 


( 11 ) 
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II, Rectangular Wave Approximation 


In order to simplify the analysis, we will assume 
a rectangular intensity distribution in the form 


I for O — r ^ a 

I = (12) 

0 for r >■ a 


Then the differential equation becomes 


d m(r) - (2trl+w) m(r ) = w o^r-a 

dr D D 


d m(r ) - w m(r) = w r >a (13) 

dr D D 


d g m(r) + _1 
dr 2 r 

d 2 m(r) + __1_ 
dr 2 r 


Letting 



2 o- I + w 


w 




D 


(Ik) 


We obtain the solution 


m 

= m^ o 4: r a 


“L 

m r > a 

c. 2 

= AI D (k 1 r) + BK (kn r) - k £ 

(15) 


*i 2 


“2 

= Cl (k.r) + X (k«r) - 1 
0 2 o 2 
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I c (x) is the modified Bessel function of the first 
kind of zero order and K q (x) is the modified Bessel func- 
tion of the second kind of zero order, A,B,C, and D are 
constants to be determined by boundary conditions, which 
can be taken as 

(1) ( a) = a) 

(2) m^( a) = ih P ( a) 

, . , x (16) 

( 3 ) 1^(0) = 0 

(ij.) ingfb) = - 1 


where b is the radius of the cylindrical absorption cell. 
Boundary condition (3) dictates B to be zero. Using 
the other boundary conditions, and equation (11), we 

L 

obtain 


2( k 2 2 ~ k l 2 ) Il(kia) 
m av = k| k l a 

I 0 (ki a ) + P ki'l^Tk^a) 
E 



( 17 ) 


T , __ I(x) is the first order modified Bessel function of the 

wnere f±rst klnd gnd 

E = k g C ^(kgaJI^Ckgb) + K l (k 2 a)I <> (k2b)] 

P= [K o (k 2 a)l 0 (k 2 b) - I 0 (k 2 a) K Q (k 2 b )3 


By definition, the saturation intensity occurs when m av 
Then we find 


1 

2 
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k l 2 - k 2 2 

kia(ki 2 -2k 2 2 ) 



It can be shown that the terms with K (k~b) are 

o 2 

negligible for reasonable values of b. For example, in 

O 

oase of SF^,, kg is about 1.17 at 1 Torr. If we take b 
to be 1 inch, then k 2 b = 2.98, However, k^ is essen- 
tially proportional to the pressure, as will be shown 
later, so that k 2 b will be smaller at smaller pressures.; 
Let us take two extreme cases: 

Case 1 . k 2 b = .1 k^a = .01 

K D (k 2 b) = 2.427 ^'(k a) =100 

^(kga) = 0 I o (k 2 b) = 1 

K 0 (k 2 a) = 50 I 0 (k 2 a) = 1 

Therefore I^Ck^a) K o (k 2 b) K^Ck^a) I Q (k 2 b) and 

K 0 (k 2 a) I Q (k 2 b) >^I o (k 2 a) K^kgb) 

Case 2. k 2 b = If. k 2 a = .!)- 

K 0 ( k 2 b ) = .011 .^(kga) = 2.18 

I x (k 2 a) = '.2 V k 2 b) = 11L * 3 

K 0 (k 2 a) = 1.115 I 0 (k 2 a) = 

Again L^kga) K (kgb) C< K^kga) I Q (k 2 b) and K Q (k 2 a) 
I 0 (k 2 b) » I Q (k 2 a) K Q (k 2 b) 

So ignoring K Q (k 2 b), we obtain 
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( 20 ) 


o 

E = I (k 2 b) K-L(k 2 a) k 2 
P = i (lcb) K (k P a) 

O d O d 

Then equation (19) can be written as 

~ k 2 2 = 1 r W + K o lk 2 a) k l 1 (£1) 

kiMk^-gka 2 ) _r r Ll i Xkia) K i (k 2 a) k 2 


Let us call the right-hand side of this equation Z, 
If we let 


(k x a) 2 _ k-i 2 

IkTaF '- 1 ? 

k 2 


( 22 ) 


¥e have 

X = 1 - 2(k^a)2 (23) 

1 — (k^a)Z 


If we plot equations (22) and (23) together, as 
a function of (k-^a) for various values of (k 2 a), the in- 
tersection point gives the value of k^a as well as X. 
Then the saturation intensity is 

I s = w [X-l] hr (2k) 


Figure 1 is the result of these cal culations. The 
intersection points determine X as a function of k 2 a. 
(X-l) versus k£ a is plotted in Figure 2. For a parti- 
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FIGURE 1 

GRAPHICAL METHOD FOR THE CALCULATION 
OF THE SATURATION INTENSITY 


9 


cular gas is equal to ,/w/D. Figure 2 shows that the 
saturation intensity changes very drastically for small 
values of k£a, For large values of k^a, the saturation 
intensity does not change much. 

We note that the theory includes the effect of the 
cylindrical absorption cell walls with the boundary con- 
dition m(b) = -1, 
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Analytical Derivation of the Absorption Coefficient 
and the Saturation Intensity 


Experimentally we normally measure the output 
power as a function of the input power, and deduce the 
variation of absorption as the intensity of the beam is 
varied. So an important quantity to know is m _ __ as a 

EL V 

function of the input power. We can write equation (17) 

6 

as 


m av = (1 - 2k ia Z) - X 
(2k x a 2) X 


(25) 


m QT . as a function of (X - 1) is plotted in Figure 

81 V 

3 for various values of k 2 &. We see that the absorption 
decreases much faster as the beam size is increased. 
Without diffusion, m 0 , r changes as 

8. V 


m 


av 


1 + I/I. 


(26) 


Differentiating this expression with respect to I 
we obtain 


m av = 1 


(27) 


Ig (1 +_!__)' 

*3 


At I = 0, we have 


™ a v lT s. - 1 


( 28 ) 
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When diffusion is included, we should obtain the 

6 


same result as I goes to zero. Now we have 


o 

m 


av ^ £ 1 " 

+ E 


uVV 1 


1 


(29) 


Approximating F/E as before and comparing equations 
(28) and (29) we obtain 


I = _w_ K(k 2 a) (30) 

2 «r ' 

where the function K (kg a) is given by 

K(k 2 a) = 1 (31) 

1 ** 2 

k 2 af I o (k 2 a) + K Q (k 2 a ) ) 

Il(k 2 a) K-j_ (k 2 a" 


w/2cr- is the multi-level saturation intensity. The func- 
tion K(k 2 a) is what we called (X - 1) before. It is the 
factor that is caused by the diffusion and wall effects, 
and is plotted in Figure 2. It is seen that the graphi- 
cal and analytical solutions give essentially the same 
results. 

Experimentally we normally measure output power 
versus input power. It is easier to do this measurement 
at high power densities. For example, if Measured = 
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FIGURE 2 

GRAPHICAL AND ANALYTICAL SOLUTIONS 
FOR THE SATURATION INTENSITY OF MOLECULAR GASES 
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FIGURE 3 

INVERSION DENSITY AS A FUNCTION OF INTENSITY 
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o o 

2.1 W/oia and l rQa ^ = 2.2 W/cm , we have $ error = 5$. 

p 

At low power densities, if l mQaSure d = .0001 W/cm and 

2 

■^real = # ®^02 W/cm., % error = $ 0% 9 However, we need to 
show how to interpret the results at high power densi- 
ties meaningfully. At large intensities, we can write 


k-, 2 = 2o- I 

D 


(32) 


This is further justified in the next section. Now 
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we can write equation (25) in the form 

d ; 

m = - u W 2 


where 


av 


r D 1 + 1*1 

= - L w 2« J 


2 cr I + 1 


W 


Z» = 2a Z 


( 33 ) 




For large values of k^, we can write 


2' = 1 


K q (k «) 


2 (k a. ) k2 


( 35 ) 


This is further justified in the next section. 
Combining equations ( 9 ), ( 32 ), and ( 35 ), we obtain 


2. 

i 


ai = 
dz 


_ K» (k 2 a) 

2 crl + 1 

w 


( 36 ) 
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where 


;« (k 0 a) = 2K 1 * k 2 a ) + 1 


k 2 a K q (k^a) 


Integrating equation (36), we find 
K. (k 2 a) « q L - In + 221 £ I oufc " hn] 


( 38 ) 


At high power densities, the first term on the 
right-hand side can be neglected. Then we have 

I out =I ln - ^L_ K'(k 2 a) « 0 L (39) 

2 cr 


Without diffusion, x-je would have 


dl = 
dz 


1 + I 


0K» 


Neglecting 1 in the denominator, we would obtain 


I , = I. - o' L I» 
out m o " A s 


(Ip.) 


Therefore, it is tempting to identify 


I ‘ = w K» k~a) 

s r 2 

2 cr 


0 ^ 2 ) 


The actual saturation intensity is given by equa- 
tion (30). Therefore if one measures I* experimentally, 

s 

the actual saturation intensity will be 
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I = [CF] I « (43) 

S 3 

where £ CF] is the correction factor. It is given by 

IGF] = K(k 2 a) (44) 

K« (k 2 a)" 

CCP] , K (k 2 a), and K' (k^a) are tabulated in Table 1. 
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TABLE 1 


Functions of Interest 



x - 1 

0.1 

75.075 

0.2 

24.39 

0.3 

13.74 

0.4 

9.191 

0.5 

6.865 

0.6 

5.473 

0.7 

4.571 

0.8 

3.935 

0.9 

3.481 

1.0 

3.132 

1.1 

2.836 

1.2 

2.639 

1.3 

2.454 

1.4 

2.325 

1.5 

2.195 

1.6 

2.092 

1.7 

2.004 

1.8 

1.928 

1.9 

1.86 

2.0 

1.802 

2.1 

i.75i 

2.2 

1.705 

2.3 

1.663 

2.4 

1.626 

2.5 

1.592 

2.6 

1.562 

2.7 

1.533 

2.8 

' 1.508 

2.9 

1.484 

3.0 

1.485 

3.1 

1.444 

3.2 

1.428 

3.3 

1.411 

3.4 

1.395 

3.5 

1.38 

3.6 

1.367 

3.7 

1.355 

3.8 

1.343 

3.9 

1.332 

4.0 

1. 321 


2 K i ( V> * i 

K^Tk^ajk^a 


82,202 
28 . 21+8 
15.845 
10.8001 
8.1672 
6.5835 
5.543 
4.811 
4.278 
3.865 
3.5323 
3.274 
3.062 
2.887 
2.729 
2. 601 
2.489 
2. 392 
2.306 
2.229 
2.16 
2.103 
2.045 
1.986 
1.949 
1.907 
1.867 
1.83 
1.795 
1.772 

1.741 
1.716 
1. 693 
1. 668 
1.647 
1. 628 
1.61 
1.59 
1.574 
1.558 


OF 


0.9133 

0,8634 

0.8672 

0.851 

0.841 

0.831 

0.825 

0.818 
0.814 
0,81 
0,803 
0.806 
0.801 
0 . 805 
0.804 

O.804 

0.805 

0.806 

0.807 

0.808 

0.811 

0.811 

0.813 

0.815 

0.817 

0.819 

0.821 

0.824 

0,827 

0.827 

0.829 

0.832 
0.833 
0.836 
0.838 
0.84 
0.842 
0.845 
0.846 
0. 848 
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IV* The Case of Inhomogeneous Broadening; 


At low pressures that are used with the saturable absorbers, 
we would expect the Doppler broadening to be dominant over 
collisional broadening. In order to see what happens in this 
case, let us rewrite equation ( 36 ) in the form 

= L( V>, k 2 a) (k?) 

Q 

ZszL + 1 

w 


where 

L(k-i a,k/3a) = 2 1 + 1 (1+6) 

k 2 a i 0 (y )k g ; k 

I^TkJa)^ K^Tk|aT 


We need to integrate equation (1+J?) in order to average 
over the classes of molecules that are broadened differently. 
The variables are oC Q , o - , and k-^. The integral is 

00 

o£. = S (tOp(tt )dr (1+7) 


h 


The letters h and i stand for homogeneous and inhomo- 
geneous respectively. If the function L(k^a, k^a) were inde- 
pendent of r 0 , then the above integral would integrate as in 
the classical case to give a result similar to equation (2). 

So let us see how this function behaves. If the beam intensity 
2 

is 1 W/cm , we would have for SFg, 
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2 crl = 1+.07 • 10^ 


at the peak of the Lorentzian gC - **©) . The relaxation rate of 

8 

SP^ at 1 mm is about 50/sec , However, we need to know the 
effective relaxation rate. It is clear, though, that 2o-I 
is much larger than w. So at intensities larger than the 
saturation intensity, it is safe to say that 

h = vpST (48) 

I© (ki a )/Ii (k-j_a) remains close to 1 for reasonable values of 

k-j_a : K^k^a) /K-j_ (k^a) is close to • 5 and does not change much 

with k a. We have 
2 

k 2 * / w (l*.9) 

k x yj 2v-I 


Letting 


= G g (*r) 


(5o) 


where C is a constant, and substituting the Lorentzian lineshape 
function g('sr)^ in equation (48^, we find 

k 2 4 / TTw f ( Av ) 2 + (-y -t ) 2 j (51) 

k l J i 

v C IAv 

g(-r) falls quite rapidly from its maximum value of 2/TTA*r, 

It is loss than one tenth of its maximum value when (y - -r ) 
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is equal to Att. With the two level model without diffusion 
we would have 

c*. = - MChr (52) 

1 + TTAy I 

gW 2 I s 

where M is the molecular density. Substituting in g(*r), we 
find 

oc h = - M Chy A -r/2TT~ (53) 

[(■r-^) 2 + ( A-r) 2 ] 4- A-r 2 I 

2 4 s 

This is broader than the lineshape function g("r) by 
J~1 + I/Ig* Assuming I/I g to be 100, the width of would 

be about 10 Ay, We see that as I increases, the width in- 
creases while the peak value decreases. 

Let us assume k^a = ,1, I Q (k^a)/^ (k^a) = 1, K Q (k^a) / 
(kg a) = .5, w = 10 2 , 2 «rl = 4 x 10 6 . Then 

/ w = 1/200 

v/ 2«r I 

Using these values, we define 

V(v)=_ 20 +1 (54) 

,01 y .25 + (-r - ~r 0 ) 2 +'.5 
Ay 2 
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We see that the denominator above changes very little 
for (•** - u* 0 ) 2 / 10. This function is plotted in 

Figure ij.. It is seen to be varying very slowly. Therefore, 
for all practical purposes, it is safe to assume that L(k^a, 
k 2 &) can be assumed to be constant while integrating to 

find Then we immediately obtain. 



Ltk^a, k 2 a) od Q 
J 1 + 2 o-l/w 


( 55 ) 


At intensities larger than w/2 o- , we can neglect the 
modified Bessel functions of the first kind, obtaining 


o = K« (k 2 a) c^ 0 
y] 1 + 20-I/W 


( 56 ) 


Integrating equation (56) gives 

.1 


K> (k 2 a)o«i 0 L = 


= In 


(1* 20-1 )9- - 1' 

w 

^ l* 2o-I )tr * 1 

w 


+ 2(l* 2q-I ) 

w 


'] 


out 


■in 


( 57 ) 


When 2vl/w >$>*1, namely at high intensities, we obtain 

K ' (k 2 a) “V 1 = ['/ T out " </hn] <58) 

In the case of the multi-level model without diffusion 
we have 
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FIGURE 4 

THE FUNCTION Y (r) 
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( 59 ) 


1 <3 1 s - °^o , 

I <V < 1 + I >2 

<S 

At high intensities, we would obtain 

^ 0 x. = _a. [ /X^ t - ] (60) 

Therefore, the saturation intensity I*,, can be defined 

■3 ’ 

as - 

II « [K‘(k 2 a)] 2 (61) 

2 o- 

The actual saturation intensity is expected to be some- 
what lower than I’ as In the homogeneous case. 

3 
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V, Summary 

The main results of the foregoing analysis can be sum- 
marized by the following equations: 

k 2 = Jw/ D (62) 

w — effective relaxation rate 
D = effective diffusion constant 

k-^ = J 2q-I + w (63) 

U~ = absorption cross-section 
I = photon intensity 
a - radius of the beam 


K(k 2 a) = 




1 - 


k 2 a f I 0 (k 2 a) + K o (k 2 a) 

L rjjk^al K-l (k^a) -» 


K* (k 2 a) = 2 K x (k 2 a) + 

k^a K q (lc 2 a) 


( 65 ) 


Ltk^a^kga) = 2 • 


k 2 a k^Ja ^o^ k 2 a) 

k ia + K 1 (k 2 a) 


+ 1 


( 66 ) 
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Homogeneous Case 


o£ = L(k x a, k 2 a) Q ^ n 


(67) 

2 o*I + 1 



w 



* = K'(k 2 a) 

I >> w 

p__ 

(68) 

2 o*I + 1 






'out ~ X in = -IJ °^o L 1 

> w/2o- 

(69) 

I' = _w_ K'(k a) 


(70) 

2<r 



i s = Lcp] i« 


(71) 

[CP] = K(k 2 a) / K»(k 2 a) 


(72) 

Inhomogeneous Case 



oi = L^a, k 2 a) c* 0 
v/ 1 + 2 o-I/w 


(73) 

= K'(k 2 a) od 0 I ^ 

w 

(74) 

^1 + 2o-l/w 

2o* 


VW “ 4 ^in ~ “ J^T^o ^/2 

1^ w 
2^ 

(75) 

I> = _w_ [K* (k 2 a)] 2 
2 


(76) 
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VI. Conclusions 


The results above indicate, first of all, that 
diffusion plays a very significant role in the 
saturation of molecular gases. Secondly, the satura- 
tion intensity in the case of inhomogeneous broaden- 
ing is seen to be quite different from the saturation 
intensity in the case of homogeneous broadening. As 
a matter of fact, using equations (69) and (75), we 
can write 


out 


- p in = - [a a ' ( k 2 a) ] A (77) 


for the homogeneous case, and 


< aK,(k 2 a » 


B (78) 

for the inhomogeneous case. A and B are constants. 

The theoretical results derived herein were 
tested with seven molecular gases using various lines 
of the CC >2 laser around 10.6 microns. The gases were 
SF 6 , BC1-, C 2 H 3 C1, NH 3 , C 2 H ]j F 2 , CH 2 CF 2 , and CCl^. 


The most interesting observation during these 
experiments was that the output power did not vary 
much for constant input power when the radius of the 
beam was changed. Going back to equations (77) and 
(78), and taking this ob- 


27 



p 

servation into consideration, we can say that if a K' (kga) 
does not vary much when a is varied, the broadening should 
be homogeneous, and if aK' (k^a) does not vary much when a 
is varied, the broadening should be inhomogeneous. Table 2 
is prepared for this purpose. It is seen that aK' (k^a) remains 
more or less constant in a large range of k 2 a whereas a K' (k^a) 
increases quite fast with k 2 a. Of course, we would expect 
aK’ (k^a) to start increasing linearly with k 2 & since K’ (k 2 a) 
approaches 1 for very large values of k 2 &. 

The only way the function of aK' (k 2 a) can be constant 
with respect to a is that K’ varies inversely with a. 

From Table 3, we find 

K« (k 2 a) = 5.15 (79) 

a 

in the range of k 2 a indicated in Table 3« This result is 
most interesting. From a physical standpoint, the variation 
of the effect of diffusion inversely with the radius of the 
beam is intuitively sensible. 

So if the inhomogeneous case is valid at low pressures, 
we can conclude that the variation of the output power with 
the radius of the beam at constant input power is negligible. 

We can write this conclusion in the form 
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TABLE 2 

TEST FOR BROADENING 


V 

k 2 a.K* (k 2 a) 

(k 2 a) 2 .K' (k 

0.1 

8,22 

O.82 

0.2 

5.65 

1.13 

0.3 

4.75 

1.43 

o.4 

4.32 

1.73 

0.5 

4.08 

2.04 

0.6 

3.95 

2.34 

0.7 

3.88 

2.72 

0.8 

3.85 

3.08 

0.9 

3.85 

3.47 

1.0 

3.87 

3.87 

1.1 

3.89 

4.23 

1.2 

3.93 

4.72 

1.3 

3.98 

5.07 

1.4 

4.04 

5.66 

1.5 

4.16 

6.24 

1.6 

4.23 

6.77 

1.7 

4.31 

7.37 

1.8 

4.38 

7.88 

1.9 

4.46 

8.4-7 

2.0 

4.54 

9.08 

2.1 

4.63 

10.72 

2.2 

4.70 

10.34 

2.3 

4.77 

H.45 

2.4 

4.87 

12.18 

2.5 

4.96 

12.90 

2.6 

5.04 

13.64 

2.7 

5.12 

14.34 

2.8 

5.21 

15.11 

2.9 

5. 32 

15.96 

3. 0 

5.4o 

16.74 

3.1 

5.49 

17.57 

3.2 

5.59 

18 .45 

3.3 

5.69 

19.27 

3.4 

5.76 

20.16 

3.5 

5.86 

21.10 

3.6 

5. 96 

22.05 

3.7 

6.04 

22.95 

3.8 

6.1 4 

23 c 95 

3.9 

6.23 

24. 92 
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(80) 


where the subscripts 1 and 2 refer to beam of radius a^> and 
the beam of radius a^. Equation (80) is expected to be valid 
at constant input power and low pressures. Using equations 
(70) and ( 76 ), we can also conclude that the saturation in- 
tensity varies inversely as the radius of the beam for the 
homogeneous broadening and varies inversely as the square of 
the radius of the beam in the case of inhomogeneous broadening, 
at low pressures and small radii. 
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VII. Appendix 


Derivation of Effective Relaxation Rates and Diffusion 
Constants 

Neglecting Diffusion, the rate equations can be written 

6 

as : 


m 2 = 


U 


u 


W 2 ”M 2 - ^_w 2 J M 2 + X % - W(M 2 - + P 2 U (A.l) 

i i 


o U 

“k 


*L " 


w k _ln k “ 2- m k + X "ikV* + p k u 


U U “I\ U 

i 1 
1 


(A e 2) 


w iS - 


”11^ + W( m 2 - 1\) + Pj 1 (A. 3) 

i i 


• 1 _ 
^k “ 


1_ 1 

™k ^k " 


+ 

i i 


“ik ”! 1 + P k 1 


(A.W 


In these equations various terms are defined as follows: 
W = stimulated transition rate between upper and lower absorp- 
tion levels. 

w^u = relaxation rate of the upper band kth level to all 
•. lower : *r band levels 
P k^ = P urn P ing rate of the kth upper level 
P, 1 = pumping rate of the kth lower level 
w^ — relaxation rate of the lower band kth level to all 
other lower band levels 

ss relaxation rate of the kth upper (lower) band lpvel 
to the ith upper (lower) band level 
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= Sg H^, where gp and are the degeneracies of levels 
2 and 1 


“k 


= population density of the upper (loiter) band kth 
level 


The number of equations and terms are actually infinite 

since k can assume any value. At steady state all the deri- 

v 

vati<yes are zero. Adding equations (A.l) and (A,2), we get 


0 






u u 

+ 



„ U u 
w ik m k 


- W(M 2 - K[) 




U.5) 


The second and third terms cancel out. Assuming all the 
pumping rates are zero except for the upper absorption level, 
we obtain 


r- X Vm" 1 * 

i 1 M ? - W(M2 - \ ) + P ? u = 0 (A. 6) 

L pu J 


Similarly, considering equations (A. 3) and A. ip) , we get 


w i lm i 1 “* 


Hi 


Mi+W(M - K^) + p 1 1 ~ 0 (A* 7) 
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We can assume all the rotational relaxation rates in 
each band are the same, equal to w u and w^, Then equations 
(A, 6) and (A. 7) are the same as the two-level model equations 
with the relaxation rate of the upper level replaced by w u 2m* u / 

i 

M 2 and the relaxation rate of the lower level replaced by 

w 1 Z.mj 1 /^. 

In order to include the effect of diffusion in the rate 
equations, we have to add DV M(c) term into equations (A,l) 
through (A.4). For example, equation (A.l) becomes 


$2 = “ w 2 Um 2 + D ^ 2m 2 



w 


u 


2i 





w 



u 

i 


- W(M2 - M[) + P 2 U 


(A. 8) 


Again assuming stead-state, if we add all the upper level 
equations, we obtain 


o 7 1*1 J ,, 

- w 2 W 2 + DV ( jr j ) H 2 - W(M2 ) + P 2 U = 0 (A.9) 

^2 


u 


So we can claim that the effective diffusion constant 
for the upper level is D^_ m^ /iig* We obtain a similar 

result for the lower level by adding the lower level equations; 

namely, the effective diffusion constant for the lower level 

1 

D, becomes D 2L- m^ /ii^, Assuming that the rotational constant 


33 



for the upper and lower bands are the same, and will 
be nearly equal. 
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